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By exploiting the concepts of magnetic group theory we show how unidirectional behavior can 
be obtained in two-dimensional magneto-photonic crystals (MOPhC) with uniform magnetization. 
This group theory approach generalizes all previous investigations of one-way MOPhCs including 
those based on the use of antiparallel magnetic domains in the elementary crystal cell. Here, the 
theoretical approach is illustrated for one MOPhC example where unidirectional behavior is obtained 
by appropriately lowering the geometrical symmetry of the elementary motifs. One-way transmission 
is numerically demonstrated for a photonic crystal slice. 

PACS numbers: 42.70.Qs, 41.20.Jb 



The unique dispersion properties of photonic crystals 
(PhC) have led over the past decades to the theoretical 
and experimental demonstration of many unusual and 
extraordinarily enhanced optical phenomena [TJ |2] . Pho- 
tonic crystals made of magneto-optical (MO) materials, 
where time reversal symmetry is broken, open a new re- 
search field on non-reciprocal behavior in optics. Several 
papers have thus proposed theoretical non-reciprocal op- 
tical demonstrators based on such MO PhCs. Indeed, 
time reversal breaking through the presence of MO ma- 
terials is greatly enhanced by exploiting the strong dis- 
persion of photonic crystals. In regions of low group 
velocity, this can lead to unidirectional freezing of light 
[3, 4 : forward light propagates in the crystal with finite 
velocity whereas backward light is stopped. Recently, Yu 
et al. have even shown that time reversal breaking can 
lead to a non-reciprocal shift of bandgap edges in a one- 
dimensional periodic MO Bragg stack [5 . In this case, 
there appears a frequency range within which light can 
freely enter the MO PhC in one sense, while it is totally 
reflected in the opposite sense. However, Figotin et al. 
have shown that optical unidirectionality requires both 
time reversal and space reversal breaking of the structure 
[4]. Till now the latter condition has only been achieved 
in theory by using periodic structures with antiparallel 
magnetic domains in the unit cell. Such structures would 
require a very precise control of magnetic domain walls. 

In this Letter we show how, by exploiting the concepts 
of magnetic group theory, one-way behavior can be ob- 
tained in a generic manner in two-dimensional MO PhCs. 
A general scheme can thus be used to classify all geome- 
tries of a magnetophotonic crystal with possible unidirec- 
tional effects. In essence unidirectionality requires first 
and foremost that the symmetry group of the periodic 
structure be reduced in such a manner that an arbitrary 



Bloch k- vector of the Brillouin zone (BZ) and its time re- 
versed opposite are not equivalent. For this, time rever- 
sal breaking is a necessary but not a sufficient condition. 
In general, one has to imagine a crystal layout in which 
none of the present symmetries allows the transformation 
k — > — k. In what follows, we illustrate how group the- 
ory concepts can be used to obtain non-reciprocity in a 
uniformly magnetized PhC by simply reducing the sym- 
metry of the elementary motifs. Theoretical predictions 
are then numerically validated by demonstrating a pro- 
nounced one-way transmission through a finite slice of 
the crystal. We will finally indicate how properly choos- 
ing a 2D MO PhC layout among all possible magnetic 
groups can generalize such a one-way behavior. 




FIG. 1: (a) Uniformly magnetized hexagonal magnetopho- 
tonic crystal with a motif of triangular symmetry. Symme- 
try elements indicated in the unit cell (dotted black lines) 
are: rotation centres of order 3 (black triangles), anti-mirror 
planes (time reversal + mirror operation - continuous red 
lines) and anti-glide planes (dashed red lines), (b) First Bril- 
louin zone (hexagon) and its irreducible part (IBZ) (shaded 
triangle Y — K — K'). Black crosses form the "star represen- 
tation" of an arbitrary k-vector. Both a Bloch mode and its 
time reversed counterpart (red cross) are within the IBZ. 

Figjlla) shows the specific (but not exclusive) layout 
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that we have chosen to illustrate the general concept. A 
clover-like airhole motif repeated on a hexagonal lattice 
is etched into a MO background material that is magne- 
tized parallel to the air holes. The possibility of achieving 
band structure non-reciprocity for this 2D MOPhC with 
uniform magnetization can be understood from the anal- 
ysis of the crystal symmetries. Because the axial vector 
M is reversed under time reversal, time reversal is not a 
symmetry operation for Maxwell's equation. As a result, 
the time-reversed Bloch mode (k — > — k) of an arbitrary 
k-vector in the Brillouin zone (BZ) will a priori not be 
a member of the k-star representation of this solution [6] 
(Figjljb)). This is a first step towards a non-reciprocal 
band structure (uo n ^ ^ tj n> _k), but other symmetry op- 
erations may still be present that transform a given k- 
vector into its opposite. In order to precisely determine 
the general k-star of the BZ, one needs to identify the 
"magnetic" space group elements of the MOPhC. The 
mathematical derivation of this space group is outside 
the scope of this paper. Using the construction derived 
by Indebom [7 one proves that the plane magnetic 2D 
group M of the MOPhC of Fig jl] is given by 

M = p3 + i(p3ml - p3) p3m'l (1) 

where t is the anti-unitary identity operator (i.e. the 
time inversion operator, (M. = — M). p3 and p3ml 
are the conventional notations for the plane triangular 
groups [8] and the prime stands for combination with 
time reversal (t). Fig[IJa) illustrates the symmetry ele- 
ments of p3m'l. The unitary rotations of order 3 (sub- 
group p3) are indicated by black triangles. The anti- 
unitary elements of the complement t(p3ml — p3) are 
combinations of time reversal and mirror operation: a 
pure mirror operation parallel to M would reverse M. 
These anti-mirror planes (and anti-glide planes) are rep- 
resented by continuous and dashed red lines respectively. 
Operating the unitary and anti-unitary point group oper- 
ators of the symmetry elements on an arbitrary k— vector 
leads to the star representation as shown in Figjljb). The 
absence of time reversal combined with the symmetry 
reduction of the elementary motif limits the number of 
branches in the k-star to 6 without inversion center in k- 
space (instead of 12 for a symmetric hexagonal lattice). 
In consequence, the band structure is non-reciprocal and 
the irreducible Brillouin zone (IBZ) is larger than that 
of a symmetric hexagonal PhC. A possible choice for the 
IBZ is the shaded area in Figjljb), which now includes 
both a general Bloch mode and its time reversed coun- 
terpart. Indeed, because of the anti-unitary character of 
the mirror symmetry elements in the p3m'l group, — k is 
equivalent to m xy \i. 

— k = tk rsj tm! xy k ~ ttm xy ~k ~ m xy h (2) 

Except for the Bloch vectors k on the Y — M line (parallel 
to the anti-mirrors in p3m'l), all k-points in the IBZ are 



intrinsically non-reciprocal. This explains the difference 
in notation between the K and K' points in the IBZ. The 
non-reciprocity and possible unidirectionality will thus 
be more pronounced the further one moves away from 
T — M. The highest effects are then expected around the 
if -point (and the non-equivalent time reversed K r point). 
Using only general arguments from magnetic group the- 
ory we have with this specific example shown how a can- 
didate MO PhC layout for one-way effects can be identi- 
fied and where in the photonic bands the strongest non- 
reciprocal effects can be found. It suffices to choose a 
suitable magnetic group that doesn't possess any inver- 
sion symmetry in k-space and then use a geometry for 
the elementary motif that is in accordance with the cho- 
sen magnetic group. As illustrated by the example, this 
can be achieved with just a slight geometric symmetry 
reduction. 

Let us now illustrate this non-reciprocity from band 
diagram calculations for the 2D MOPhC of Figjl] The 
permittivity tensor of the MO material, for instance Bis- 
muth iron garnet (BIG, BisFesO]^) is written as 

/ e +jg z -jg y \ 
e = e -jg z e +jg x (3) 
\+39y ~J9x e / 

where e is the isotropic permittivity of the material, 
with g xi g y andg z the components of the gyrotropy vector 
g = #1m, with 1m the unit vector along the magneti- 
zation direction, and g is the gyrotropy constant. This 
tensor form is an approximation that leaves out the crys- 
talline description of the material, and the non-reciprocal 
anisotropy only stems from the direction of the magne- 
tization. For most MO materials, this is an acceptable 
description. The gyrotropy constant g is related to the 
specific Faraday rotation through F = We take 
e = 5.0625 and g = 0.1 for BIG at near infrared fre- 
quencies [9 . The antisymmetry of the off-diagonal ele- 
ments, i.e. the source of the non-reciprocity, is imposed 
by Onsager's principle (e/^(M) = M)), while power 
conservation for transparent materials (e and g G R) im- 
poses the hermiticity of e. The sign of the gyrotropy 
vector components g x , g y and g z changes under time re- 
versal since they are proportional to the magnetization 
unit vector. Besides, any mirror plane parallel to z will 
induce a change of sign of the g z . For the 2D crystal of 
Fig|TJ g x = g y = 0, and g z = g is a constant throughout 
the uniformly magnetized MO material. Within these 
conditions, the 2D MOPhC is said to be in Voigt config- 
uration [10] . The specific form of the permittivity tensor 
leaves the TE (E z = 0) and TM (E x = E y = 0) Bloch 
modes of the 2D crystal uncoupled, while it leads to a 
non-reciprocal correction of the TE mode dispersion via 
the antisymmetric elements (e xy = — e yx ). 

Fig. [2|a) shows isofrequency contours in the recipro- 
cal space of the crystal for normalized frequencies of the 
second TE band near the K and K r points: 0.34 ^ u*(= 



3 



g = +0.5 g=0 g=-0.5 




k x (x a/27r) k x (x a/27c) K (x a/2;c) 



FIG. 2: (a) TE isofrequency contours for 0.43 ^ to* ^ 0.438 in reciprocal space for the MOPhC of Figjlfa) with, from left to 
right, g = 0.5, 0, —0.5 . The strong non-reciprocity of the band structure around the K(K') points induces a local forbidden 
zone for Bloch vectors with k y ~ |j for g > or the opposite for g < 0. For a nonmagnetic PhC (g = 0) with reduced-symmetry 
motif, the band structure preserves the hexagonal symmetry of the lattice and is reciprocal, (b) Principle of a unidirectional 
uniformly magnetized MOPhC mirror based on the one way k— space band gap. A TE plane wave (E z = 0) with k y « |^ 
and ujk ^ u>* ^ wk 1 (second TE band) will be totally reflected in one sense and partially transmitted in the opposite sense. 
The condition for the optimum incidence angle depends both on the frequency and the index of the incidence medium via 
uj*nin C sin (ft) = 1/3. 



ooa/2-Kc) ^ 0.348. These contours have been obtained us- 
ing a 128 x 128-planewave expansion of the unit cell of the 
periodic permittivity tensor profile pJU . The geometrical 
parameters of the crystal are r/a = 0.2 and d/a = 0.2. 
The three graphs differ only by the sign and/or value of 
the gyrotropy coefficient g. The middle graph shows the 
contours when the external magnetic field is switched off 
(g = 0). As expected, the bands are reciprocal and the 
isofrequency contours exhibit the hexagonal symmetry of 
the Bravais lattice, notwithstanding the reduced triangu- 
lar symmetry of the motifs. On the contrary, as soon as 
the magnetization is switched on (left and right graphs 
in Fig. [2|a)), the time-reversal symmetry is broken, and 
the reduced symmetry of the motif manifests itself clearly 
in the shape of the isofrequency contours. Incidentally, 
the left and the right graph, obtained for exactly oppo- 
site values of g, transform into each other by a /c-space 
inversion, which is equivalent to time reversal. This cor- 
roborates the validity of our numerical method. Note 
that to highlight the non-reciprocity of the bands we took 
g = ±0.5, which is about five times greater than a realis- 
tic material constant. The most striking consequence of 
the time-reversal symmetry breaking is the relative fre- 
quency shift of the second TE band near the K and K' 
points. For g > 0, the frequency at K is higher than the 
one at K the situation is reversed for g < 0. This leads 
to the formation of local forbidden zones in the recipro- 
cal space: light propagation is forbidden at normalized 
frequencies close to ou* = 0.4385 and wavevector projec- 



tions on T — K close to k y = |^ (if g > 0) and — |j (if 
g<0). 

A one-way mirror operation can thus be obtained for 
an incident TE plane wave at these frequencies provided 
that the interface between the photonic crystal and the 
incidence medium is aligned in the Y — K direction and 
that the k y component of the incident plane wave falls 
within the forbidden zone. The same plane wave incident 
from the opposite side under the same angle (i.e. with 
opposite k y ) will have finite transmission. This princi- 
ple is schematized in Figj2^b). The angular acceptance 
of the mirror is directly related to the width Ak y of the 
unidirectional gap in k— space. This width is itself de- 
pendent on the operation frequency and reaches its max- 
imum at the nonreciprocal midgap frequency (halfway 
between the frequencies at the K resp. K' points). The 
useful bandwidth of the device can be identified with 
this maximum width. For realistic material parame- 
ters (g = 0.1), one obtains Ak* = Ak y a/2ir « 10" 3 at 
the midgap frequency. The normalized frequency band- 
width Auu* is of the order of 10 -4 . In a pessimistic 
approximation, the MO PhC slice thickness required 
for good non-reciprocal extinction is estimated to be: 
L ~ 1/Aky w 150a, which represents about 40A at 
MO gap frequencies. The device bandwidth for a non- 
reciprocal extinction ratio larger than 100 is estimated 
to be about 100GHz at near infrared frequencies. The 
unidirectional transmission of the mirror was verified us- 
ing an in-house developed 2D FDTD (Finite Difference 
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FIG. 3: (a): Normalized power transmission spectrum of the 40 period thick uniformly magnetized unidirectional MOPhC 
mirror calculated by FDTD in the frequency range of interest. The gray curve shows the non-reciprocal extinction ratio. A 
relative bandwidth « 10 -3 is obtained for a one-way extinction ratio larger than 40. (b): H z steady-state field maps 
under TE plane wave illumination in forward and backward direction at a;* = 0.457. In the forward direction light is resonantly 
transmitted while in the backward direction the unidirectional band gap rejects the incident light. 



Time Domain Method) model with a tensorial descrip- 
tion of the dielectric permittivity [12]. The results of the 
FDTD calculations are shown in Figj3] The mirror con- 
sists of a finite 40-period thick MOPhC slab embedded 
in a homogeneous isotropic medium with e equal to that 
of the BIG background of the MOPhC. Both the PhC 
and MO material parameters are the same as in Fig. [2] 
except for g whose value (—0.1) is more realistic. The 
mirror is tilted at an angle of 20° (in accordance with the 
condition of Figj2^b)), and is illuminated by a TE plane 
wave. Figj3^a) shows the normalized power transmission 
spectrum under forward and backward illumination and 
with a uniform — z magnetization (g = —0.1). This cor- 
responds to the situation in the right graph of Figj2|a). 
The non-reciprocal transmission is clearly observed. A 
one-way transmission with a non-reciprocal power ex- 
tinction ratio close to 400 is obtained in the presence 
of uniform magnetization. It may be noticed that the 
normalized midgap frequency is about 2% higher than 
that calculated from the plane wave expansion method 
(Fig{2|. This discrepancy is believed to be due to dif- 
ferent resolutions used in the two methods. A visual 
illustration of the unidirectional light rejection is given 
in Figj3^b-c), which shows the H z field distribution of 
a continuous (hj* = 0.457) TE plane wave incident ei- 
ther in the forward or the backward direction. The light 
rejection in the backward direction is striking. 

In conclusion, we have demonstrated a general ap- 
proach to obtain 2D MO PhC structures that present 
optical one-way gaps. Using concepts of magnetic group 
theory we have thus proven that strong non-reciprocal 
on-off properties are not necessarily limited to structures 
with antiparallel magnetizations in the unit cell. Instead 
of that, the use of a motif with slightly reduced symme- 
try suffices to achieve strong nonreciprocal light rejection 
in the presence of uniform magnetization. The use of an- 
tiparallel magnetic domains is just another specific case 
among the possible realizations of a non-reciprocal mag- 



netic group. Shubnikov has shown that in two dimensions 
magnetic periodic structures can be classified into 46 dis- 
tinct groups [13]. Among these only five groups do not 
allow an inversion of k-space with uniform magnetiza- 
tion. One of these five groups is the p3m f l investigated 
in this work. From a fundamental point of view this work 
is thus believed to open the way to a generic scheme for 
one-way optics in 2D photonic crystals. 
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